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■ Abstract 



A theoretical study of single electron capacitance spectroscopy in quantum 



dots is presented. Exact diagonalizations and the unrestricted Hartree-Fock 

o ' 

approximation have been used to shed light over some of the unresolved as- 
pects. The addition spectra of up to 15 electrons is obtained and compared 
with the experiment. We show evidence for understanding the decrease of the 
single electron tunneling rates in terms of the behavior of the uj — > spectral 
weight function. 

Single electron capacitance spectroscopy (SECS) Jl],|2| has been a breakthrough in the 
experimental knowledge of the electronic structure of a quantum dot (QD). Ashoori and 
co-workers |1]|J have been able to determine the energies required to introduce electrons one 
by one, from to 50, into a QD. The electrons tunnel into the QD by means of a vertical gate 
bias, and change the capacitance of the device. The measurement of that capacitance as a 
function of the Fermi energy Ep in one electrode shows a discrete set of almost equally spaced 
peaks of different intensities. A peak appears whenever Ep = fi(N) = Eq(N) — Eq(N — 1), 
with Eq(N) being the ground state (GS) energy of iV electrons in the QD. In this way, 



they were able to measure not only the addition spectra as a function of N, but also the 
tunneling rates of single electrons entering the QD when the energy balance is fulfilled. The 
main output of the experiment is to reveal a very rich dependence of /J*(N) on an applied 
magnetic field B as well as to show strong suppressions of the tunneling rates in well-defined 
regions of the B — N space. The constant interaction (CI) model [[§ for the electron-electron 
interaction cannot explain many of the features observed in /i(iV) vs. B, the more serious 
drawback of the CI model being that a single-particle picture cannot account for the observed 
suppression of the tunneling rates in this experimental set up. 

Theoretical attempts to go beyond the CI model and to explain the properties of elec- 
trons in QD's run from the simplest Thomas- Fermi ||, Hartree 0] and Hartree-Fock || 
descriptions of the electron-electron interaction to the more sophisticated exact diagonal- 
izations of the complete Hamiltonian p|-|TT|1. So far, these "exact" calculations have been 
performed in two regimes: either for small number of electrons (2 or 3), or restricted to 
the Hilbert subspace of the two first Landau levels ,i.e., first Landau index, spins up and 
down. The aim of this work is the understanding of the SECS experimental results which 
requires two ingredients lacking in previous theoretical attemps: i) the analysis of fi and 
wavefunctions in a broad range of B and iV and ii) the study of peak intensities controlled 
by tunneling rates. 

The Hamiltonian for interacting electrons is 

H = £ e i,A,a C i,o + j 12 V ijk icl a c] al C k ^Ci )(7 (1) 
i,a i,j,k,l,a,a' 

where the sum runs over the single-particle states of a two-dimensional parabolic poten- 
tial defined in the xy plane, characterized by a frequency Uq, and embedded in a magnetic 
field, perpendicular to the xy plane. ej jCr is the single-particle energy of a state i with spin a 
(an infinitesimal Zeeman term is included to remove degeneracies in spin multiplets). These 
states are characterized by the Landau index n (0 — > oo), and by the third component of 



the angular momentum, m (—n — > oo) ||11|| . Vijki are the Coulomb interaction terms given 



by four-dimensional integrals which include a cut-off (A c ) to account for the finite extension 



of the wave functions in the z direction, and and c are the creation and annihilation 
operators. Throughout this work, the standard values of the effective mass and dielectric 
constant in GaAs have been taken. u is given by huj = 5.4 meV and A c = 10A (this choice 
is discussed below). 

An exact diagonalization of the Hamiltonian up to 5 electrons is undertaken. Due to the 
circular symmetry, the diagonalization can be done in separate subspaces of configurations 
with the same z-component of the total angular momentum, M, and ^-component of the 
total spin, S z . The lowest-lying eigenvalues and eigenfunctions in all the different subspaces 
were obtained using standard diagonalization routines and Lanczos-like procedures, and the 
GS is picked among all of them. A relative error less than 0.1 % for the energy in all the 
cases was routinely achieved by including enough n and m indices of the single-particle basis 
set: Up to n = 1 for high fields and n = 2 below 1 Tesla. The evolution of fi(N) with B 
for up to iV = 5 is shown in Fig. [I] (solid lines). For comparison, we show the results 
obtained with a basis set restricted to the n = Landau index (dashed line) |[L1|| . This case 
is much less time consuming than the previous one, but significant deviations in the n(N) 
plot appear even at large values of B, as N increases. From now on, the discussion is based 
on the more-than-l-Landau-index results. 

Solid dots in Fig. |I| mark the points at which the complete polarization of the electrons 
in the QD occurs. A small bump in n(N) vs. B is observed at those points, but at values 
of B higher than those found in the experiment of Ref . || . There is an overall disagreement 
between the scales of Fig. [l| and those of the data obtained in Ref. 0. The confinement 
due to the external potential and characterized by ujq sets the curve /i(l) vs. B. Once this 
is fixed, the electron-electron interaction determines the separation /x(2) — at B = 0. 
Our calculated value is a factor of almost 2 too large compared with the experimental one. 
This might indicate that the value of A c chosen in the calculation of V^u is too small, giving 
a too strong interaction. On the other hand, the bump in /x(2) indicating the singlet-triplet 
transition lies at a higher magnetic field than in the experimental data. One should increase 
the interaction, i.e. shorten the cut-off, in order for this transition to take place at lower 
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fields (a realistic Zeeman term does not affect this transition noticeably). The only way of 
improving both scales at the same time is to weaken the strength of the interaction and to 
lower the value of ujq\ It is the dependence of on B that will be now in disagreement 
with the experiment. All this seems to indicate that other effects like non-parabolicity or 



disorder may be important [yj|. For simplicity's sake, we have chosen to work with the 
parabolic confinement that gives us in agreement with the experiment. Bearing in 
mind the scale disagreements between theory and experiment, this choice does not affect 
the main conclusions of our work. 

For what follows, it is convenient to define a "bulk filling factor" v in a QD as v = 
J2 n a v n-n,<j-, with v n<m ^ being the single-particle state occupations. This "bulk filling factor" 
reflects the Landau levels occupation close to the center of the QD as B — > oo, and it is 
consistent with the usual definition of v for a two-dimensional system. With this convention 
for v, it is possible to identify a phase diagram on Fig. [I]. The region to the left of the open 
dots corresponds to v — 2, i.e., the electrons form a paramagnetic compact droplet in the 
center of the QD, having \S Z \ minimum. Defining v m = J2 n ,a v n,m,a-, this region corresponds, 
e.g., for N = 4, to u m 2, v\ ~ 2, v 2 ~ 0, etc. A ferromagnetic state of v = 1 occurs in the 
region to the right of the solid dots, with the electrons forming a compact droplet of \S Z \ 
maximum (for N = 4: Vq ~ 1, Vx « 1, v 2 ~ 1, ^3 ~ 1, ^4 ~ ...). At even higher magnetic 
fields, a v < 1 region appears |L1], but it is not the aim of this paper to cover such a regime. 



In the region delimited by the solid and open dots, the transition from v — 2 to v — 1 
or, in other words, the paramagnetic-ferromagnetic transition (a process already considered 
in references and fLlp, takes place. First of all, the GS "melts", i.e., the charge spreads 
over a large number of single-particle states until a new compact state of total spin \S Z \ + 1 is 
formed (one can say that a spin-flip process for an electron has taken place). On increasing 
further the value of B, the same process occurs for the remaining spin-down electrons until 
complete polarization is achieved. The last spin-flip process for iV = 5 is depicted in Figure 
H]. It is in this region that the lowest-energy many-body states of different subspaces are 
very close in energy and cross each other continuously as B changes. The GS values of 



M and S z change very quickly with B in the following way: M always increases with B, 
but | S z | has usually its minimum possible value when the GS becomes one of the "melted" 
ones. On the other hand, the GS's in regions with v — 2 and v = 1 are very stable in B 
due to the large separation in energy between the GS and the rest of lowest-lying states 
of different subspaces. One can say that a "compact" QD develops in these regions. The 
"melting" process here discussed is a typical correlation effect that cannot be described 
by simple approaches as CI, Hartree or Hartree-Fock in which the parameters v m can only 
assume integer values. Although a Thomas-Fermi approach does not have, in principle, such 
a restriction, it has been used || assuming that the charge in each Landau level is quantized, 
something far from true in the "melted" regime. 

As mentioned in the introduction, not only the fulfilment of the energy balance deter- 
mines the tunneling process. For temperatures and biases smaller than the excitation energy 
(typically higher than 0.1 meV in our case) in the QD, only GS's are essentially involved in 
the process, and tunneling rates for single electrons are proportional to the spectral weight 
at uj — > 0, i.e., E F — > fi + (N) fl3f , given by 

A(iV) = E(< ) l4l<- 1) )(<- 1) |Q,.|< ) ) (2) 

i,cr 

with representing the iV-electron many-body GS wave function. Figure ^| shows 
A(iV) for iV = 2, 3, 4 and 5 as a function of B. When the electrons form a compact system 
of v — 2, A(N) reaches a value close to 1. While the polarization process is taking place in 
the QD either for N or N — 1, A(iV) is significantly reduced as a consequence of the strong 
correlation in any of the N or iV — 1 GS's. Sometimes A(iV) can even be zero due to the 
spin selection rule \S Z (N) — S Z (N — 1)| > 1/2 (in our calculation for B ~ 4.8T for N = 4 
and B ~ 4.2 and 5.4T for N — 5), It must be pointed out that, as iV increases, the interval 
of decreased A(iV) widens [see Fig. |3|, for iV = 3,4,5]. This is due to the fact that the 
polarization requires a longer variation of B and, therefore, the presence of "melted" GS's 
becomes more probable for high N. As B increases further and the ferromagnetic state of 
v — 1 is reached for both iV and N — 1, the value of the spectral weight again rises up to 1. 



This phenomenon has also been observed in the SECS experiment 0, although for a large 
number of electrons in the QD: As B increases and the QD GS turns from the "compact" 
state of v = 2 into that of v = 1, there is a region in which the capacitance signal reduces 
due to the decrease of the tunneling rate. 

A similar bump in n(N) vs. B and reduction of signal have been experimentally observed 
at lower values of B and for a large N ||. The transition from v — 2 to v — 1 could be 
studied by exact diagonalization, but this is not the case any more for a large number 
of particles. In order to go beyond the small- N regime, we use an unrestricted Hartree- 
Fock (UHF) approximation. The UHF approach is undertaken by substitution of pairs 
of electron operators by their mean-field magnitudes (c\ a ,Cj tU ), which must be calculated 
self-consistently [[14[ . 



We have performed calculation similar to the one presented in Fig. [I], with the same 
basis set used there, in the UHF approach. The UHF GS's present a similar evolution with 
B than those obtained in the exact calculation. At low field, a v = 2 paramagnetic UHF GS 
appears. In this case, the relative error of the UHF GS energy to the exact GS energy is 6 %, 
4 %, 3.2 %, and 1.8 % for = 2, 3, 4, and 5, respectively. At high field, a v = 1 ferromagnetic 
UHF GS develops. The relative error is now 0.5 %, 0.4 %, 0.3 %, and 0.2 % for N = 2,3, 4, 
and 5, respectively. In this approach, the transition from v = 2 to v = 1 takes place by 
sudden spin flips of the spin-down electrons from the center to the edge of the QD. The UHF 
approach does always give integer values for v m and, therefore, it cannot give the "melted" 
states which appear in the exact calculations. Although the UHF approximation describes 
poorly the many-body GS wave functions in the v = 2 — > v = 1 transition, the relative 
error in the GS energy remains below 6 % even in the worst cases, and tends to decrease 
with increasing N. Figure ^ shows the results of the UHF fi(N) for up to 15 electrons with 
a basis set conveniently increased up to 64 single-particle states. An approximate phase 
diagram can be traced on this plot: (a) Solids dots follow the ferromagnetic transition, i.e., 
the coming up of the v = 1 QD (as can be seen in Fig. |J it comes closer to that of the exact 
results as increases), (b) open dots, the beginning of the polarization, i.e. the beginning 



of the transition from v — 2 to v — 1, and (c) solid triangles, the appearance of the v = 2 
QD (this last transition appears also for /x(5) in the exact calculation). As the number of 
electrons increases, a bump develops more clearly (Fig. [|) and can be associated with the 
bump observed in the SECS experiment [[| at low fields and high number of particles. In 
the transition from v > 2 to v = 2 there is a new reduction of the capacitance signal in the 
experiment. This fact can be understood if a decrease of the u — > spectral weight takes 
place (as for the v — 2 — > v — 1 transition), but, unfortunately, the spectral weight is not 
well described within the UHF approach if "melted" states are due to appear. However, 
UHF describes properly the following experimental features: (a) The monotonic behavior 
of n{N) after the bump (due to the presence of the stable "compact" v = 2 state), (b) the 
oscillating behavior in the v = 4 to v = 2 transition (due to the progressive depopulation 
of higher Landau levels), and (c) the anticrossings between successive traces of n{N) in this 
region. Some of these points are given by simple models 0] while others are not, for instance 
point (c) cannot be explained by a CI model. 

A last comment remains to be done. If a finite Zeeman term is included, the correlated 
GS's that appear in the polarization process of the QD and carrying a small value of \S Z \, 
tend to be replaced by "compact" states with higher values of | *S r ^ | (as in the UHF approach) 
[ 11 1 . However, with a realistic g-Lande factor, and for a perpendicular magnetic field, the 
results presented in this work do not change noticeably. 

In this letter we have done a detailed study of the recent SECS experimental data in QD's 
@. The different regions appearing in the addition spectra as the magnetic field changes 
have been obtained and analyzed by means of exact diagonalizations and an UHF approach. 
The decrease of the tunneling rates observed in the experiment is discussed in terms of the 
reduction that the exact calculations show for the spectral weight for uj — > when the QD 
becomes "melted". 
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FIGURES 

FIG. 1. n(N) vs. magnetic field up to 5 electrons. The dashed lines correspond to the n = 
approach and the solid lines to the exact calculations. The phase diagram delimited by the dots is 
discussed in the text. 

FIG. 2. Single-particle occupations of the GS for N = 5 and for three different values of B: (a) 
Before the last spin-down electron leaves the center of the dot (B = 5 T), (b) when the transit is 
taking place through a "melted" state (B = 5.5 T), and (c) when the final spin-polarized GS is 
reached (B = 6 T). As can be seen from the single-particle occupations in (b), the intermediate 
"melted" state needs a high participation of different configurations to be described. 



FIG. 3. to — > spectral weight for iV = 2, 3, 4 and 5 as a function of the magnetic field. 

FIG. 4. n(N) vs. magnetic field up to 15 electrons in the UHF approach. The phase diagram 
delimited by the symbols is discussed in the text. The exact results for fi(N) are plotted in dotted 
lines as well as the paramagnetic- ferromagnetic transition for this case (upward-pointing arrows). 
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